Electronic excitations and the tunneling spectra of metallic nanograins 
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Tunneling-induced electronic excitations in a metallic nanograin are classified in terms of gener- 
ations: subspaces of excitations containing a specific number of electron-hole pairs. This yields a 
hierarchy of populated excited states of the nanograin that strongly depends on (a) the available 
electronic energy levels; and (b) the ratio between the electronic relaxation rate within the nano- 
grain and the bottleneck rate for tunneling transitions. To study the response of the electronic 
energy level structure of the nanograin to the excitations, and its signature in the tunneling spec- 
trum, we propose a microscopic mean-field theory. Two main features emerge when considering an 
Al nanograin coated with Al oxide: (i) The electronic energy response fluctuates strongly in the 
presence of disorder, from level to level and excitation to excitation. Such fluctuations produce a 
dramatic sample dependence of the tunneling spectra. On the other hand, for excitations that are 
energetically accessible at low applied bias voltages, the magnitude of the response, reflected in the 
renormalization of the single-electron energy levels, is smaller than the average spacing between 
energy levels, (ii) If the tunneling and electronic relaxation time scales are such as to admit a 
significant non-equilibrium population of the excited nanoparticle states, it should be possible to 
realize much higher spectral densities of resonances than have been observed to date in such devices. 
These resonances arise from tunneling into ground-state and excited electronic energy levels, as well 
as from charge fluctuations present during tunneling. 

PACS numbers: 73.22.-f, 73.22. Dj 
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I. INTRODUCTION 

In the late '90s it was demonstrated that single quan- 
tum level tunneling spectroscopy is a powerful tool for 
studying the physics of simple, noble, and magnetic nano- 
scale metals. Ralph, Black and Tinkham first used this 
technique to study the electronic energy level structure 
of individual oxide-coated aluminum nanoparticles^, and 
it was soon applied to Co, Au, Ag, and Cu nanograins 
as weli2i2i4. Even after several theoretical studies^Si&S, 
a key aspect of the data (the unexpectedly high den- 
sity of resonances in the tunneling spectra of the metal 
nanograins) remains not fully understood, as was dis- 
cussed in Ref. [2. More intriguingly, these resonances 
appear in clusters. As the nanograins are expected to 
present strong surface disorder, this bunching of reso- 
nances seemed to collide with the predictions of ran- 
dom matrix theory for the single-particle energy levels 
of small disordered conductors^. In order to explain 
these experimental features, Agam and co-workers^ ar- 
gued that the electron tunneling may occur under condi- 
tions far from equilibrium which would result in the pres- 
ence of a large number of resonances and bunching. They 
offered a phenomenological model able to account for 
some of the observed features. More recently, however, 
Davidovic and Tinkham presented alternate scenarios — 
charge traping, single occupancy of Kramers doublets, 
and non-equilibrium effects — that might be responsible 
for the bunching of the tunneling resonances that they 
observed in Au nanoparticlcs 2 . Furthermore, the present 
authors developed a microscopic model for the tunnel- 
ing spectroscopy of these nanoparticles and argued that 
clustering of resonances should arise naturally in metal 



nanograms^, even in the absence of the non-equilibrium 
transport effects introduced by Agam et alJi It has also 
been shown recently that charge-fluctuations that are 
present during the electron tunneling should generate ad- 
ditional tunneling resonances^ not considered in previous 
theories, and that penetration of environmental electric 
fields into metal nanoparticles can have striking effects 
on their tunneling spectral 

In this paper, we classify the electronic excitations that 
may take place during electron tunneling within an ultra- 
small metallic grain in terms of generations: subspaces of 
excitations with a specific number of electron-hole pairs. 
Furthermore, we propose a general microscopic mean- 
field model to calculate the quasi-particle energy levels 
in the presence of an excitation; such a model has not 
been previously presented in the literature, to the best 
of our knowledge. The generations form a hierarchy that 
strongly depends on (i) the number of available electronic 
energy levels for tunneling, and (ii) the ratio between 
the electronic relaxation rate within the nanograin and 
the bottleneck rate for tunneling transitions. The lat- 
ter is quite sensitive to the thickness of the tunneling 
barriers. The applied bias voltage, and characteristics 
of the device — grain-lead capacitances, charging energy, 
and typical energy spacing between electronic energy lev- 
els in the nanograin — determine the number of electronic 
configurations that are present in a generation. Finally, 
we present detailed results of the response of the elec- 
tronic energy structure to the excitations in an ultra- 
small aluminum grain coated with Al oxide. We find that: 
(i) The renormalization of the electronic energy levels in 
the ground states fluctuates strongly from level to level, 
and excitation to excitation due to the stochasticity of the 
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confined electronic wavefunctions that is imposed by the 
disorder present in the grain; these fluctuations depend 
dramatically on the specific realization of the disorder, 
(ii) The average single-electron energy level spacing in the 
excited states remains nearly unchanged when compared 
with the value obtained in the ground-state; changes are 
within a few percent, (iii) If the kinetics of the tunneling 
transitions, and intragrain electronic relaxation, are such 
as to admit a significant population of excited states, 
then the predicted number of tunneling resonances can 
be much higher than has been observed to date in tun- 
neling experiments in non-magnetic nanoparticles. (iv) 
Disorder is responsible for a strong sample-dependence of 
(a) the number of resonances present, as clusters, in the 
tunneling spectrum, and (b) the spectral width of such 
clusters as a function of the applied bias voltage. 



II. MECHANISM FOR ELECTRONIC 
EXCITATIONS AND THEIR SIGNATURE IN 
NONEQUILIBRIUM TUNNELING 
SPECTROSCOPY 

We are interested in the electron tunneling regime 
where: (i) the Coulomb charging energy 11 (U) of the 
grain that results from adding or removing a valence elec- 
tron, is greater than the expected average particle-in-a- 
box level spacing (5) around the Fermi energy. This is 
often the experimental situatior.I*2i&4. (ii) The electronic 
relaxation rate (r r ) within the grain is much smaller 
than the bottleneck tunneling rate (r), and (iii) tun- 
neling leaves the grain with a surplus or deficit of one 
electron with respect to the no valence electrons that 
make the grain neutral. In this regime, after an elec- 
tron tunnels in and out, the grain may be left in an ex- 
cited electronic state while still being neutral. There- 
fore, charge transport may take place via excited (non- 
equilibrium) nanoparticle states. This mechanism for 
generating electronic excitations, that may be reflected 
in the energy spectra of metallic nanograins, was first 
exploited by Agam, Wingreen, Altshuler, Ralph and 
Tinkham 5 (AWARTi) to explain tunneling spectroscopy 
experiments in Al nanoparticles. However, the original 
AWARTi model was spinless and, more importantly, the 
effects of the electronic excitations on the electronic en- 
ergy levels were treated only within the framework of 
random matrix theory. 

Let us now discuss the genesis of the tunneling-induced 
electronic excitations. We begin by considering a two- 
terminal device consisting of macroscopic source (S) and 
drain (D) electrodes separated by a thin insulating layer 
from a metallic nanoparticle (d). Following the Or- 
thodox theory of Coulomb blockade^ the electrochem- 
ical potential of the source (drain) electrode is set to: 
(j, s W(V) =E F + (-)(C D(s) /C s )eV; where V is the ap- 
plied bias voltage, e the magnitude of the electron charge, 
and Ce = Cs + Co with Cs(d) the capacitance between 
electrode S(D) and the nanoparticle. The quantum na- 



ture of the nanoparticle is taken into account through its 
discrete electronic structure. For definiteness, we shall 
assume the nanoparticle to have spin-degenerate single- 
electron energy levels \ip a ) with energy E a , and a fully 
(doubly) occupied Fermi level (\iPf)) a t energy Ep [see 
Fig. ^a)]. 12 Furthermore, by adopting Cd/Cs > 1 the 
onset of tunneling corresponds to the injection of an elec- 
tron from S into \ipF+i), with a threshold bias voltage 
given by V l s t d = (l/e){C^/C D )(U + E F+1 - E F ). As 
the number of electrons inside the grain increases by one, 
the single-electron energy levels of the grain are renor- 
malized upwards by U, within the constant interaction 
approximation^ [see Fig. ^a)]. Therefore, the subse- 
quent tunneling processes that can take place are: (i) 
return of the additional electron from level \?Pf+i} to 
S or its transmission to D; or (ii) tunneling to D of 
one of the electrons populating the single-particle lev- 
els whose energy satisfies E a + U > fi D {V^ d ). In the 
latter case, once the electron is ejected from the grain, 
the remaining electrons arc left in an excited state (|Ai)) 
that corresponds to creating one electron-hole (eh) pair 
on the ground-state electronic configuration (|G)) of the 
nanoparticle; as exemplified in Fig. Qla). Note that there 
are different possible leh-pair electronic excitations de- 
pending on which electron is ejected out. (We suppose 
that no shake-upM occurs due to the creation of the leh- 
pair excitation.) The slow electronic relaxation within 
the grain, compared with the bottleneck tunneling rate 
(recall the assumption that r/r r ^ 1) makes it possible 
to generate another family of excitations, namely, 2eh- 
pair excitations [see Fig. IHb)]. The latter, however, re- 
quire two consecutive tunneling events for their creation. 
It should be noted that the number of possible 2eh-pair 
states is greater than that of leh-pair states. 

It is important to notice that the multi-eh-pair excita- 
tions only arise as a consequence of multiple sequential 
tunneling, as described above. Hence, excited states are 
likely to exist only if T r is many times smaller than P. 
Within this context, it is convenient to introduce the con- 
cept of generations of excitations: We define a generation 
to be the set of electronic configurations that contain a 
specific number of eh-pair excitations due to consecutive 
electron tunneling in and out of the grain. Hence, gen- 
eration n (Gjj=n) contains the subset of n eh-pair excita- 
tions that arise after n such consecutive pairs of tunneling 
events. Figure ^b) shows examples of members of gen- 
eration 1 and generation 2, indicated by G#l and G#2. 
It should be clear that the generations constitute a hier- 
achy: It is not possible to have an element of generation- n 
without having generated previously elements of its n— 1, 
n— 2, n — 3, • • • , 1 ancestors. Hence, in order for n eh-pair 
excitations to occur with significant probability the bot- 
tleneck tunneling time (r _1 ) should be nearly 2n times 
smaller than the electronic relaxation time (r^ 1 ). Fi- 
nally, the applied bias voltage determines the maximum 
number of generations that may be present in a device, as 
it constrains the number of available energy levels: The 
maximum possible number of eh pairs is 2n if n levels are 
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available in the grain for tunneling; therefore, at most 2n 
generations might exist. 

To find the total number of excitations present in gen- 
eration n we note that, contrary to the case in closed 
electronic systems, the tunneling-induced electron-hole 
pair excitations may have a net total spin polarization 
as large as S z — S e z + S% ; where St^ is the largest 
possible total spin polarization of the excited electrons 
(holes). This leads to a high multiplicity of the elec- 
tronic configurations that may be present in generation 
n, as these configurations result from exhausting all the 
possible electron and hole arrangements once the restric- 
tion S z = is waived. If N (M) orbitals are available 
to accommodate the excited electrons (holes), the total 
number of electronic configurations up to generation n 

is: E"=o {I ){1 M )\ where 2n < N and the factor of 2 
accounts for spin polarization. 

The above mechanism for generating electronic exci- 
tations within a metallic grain in the Coulomb blockade 
regime, and the subsequent identification of a hierarchy 
for these excitations constitutes a generalization of the 
AWARTi model. In order to invoke this mechanism as 
a possible interpretation of experimental tunneling spec- 
troscopy data, one has to realize that the one-electron 
levels in the ground-state of the grain are not spectrally 
rigid. - 14 In particular, a net electronic interaction due to 
the redistribution of the charge within the neutral grain 
due to the excitation renormalizes the one-electron lev- 
els: The ground-state single-electron energy level struc- 
ture (\ip a ),E a ) changes to (\ip a )x n , E aX )j with \ X n) an 
excited state in generation n. Figure Qjc) schematically 
shows this renormalization in a particular case (see Fig. [21 
also). In general, for an excitation with a given number of 
eh pairs, the renormalization of level \ip a ) depends on the 
the specifics of the excited state, and therefore exhibits 
fluctuations] see Fig. [2] This is due to the stochastic- 
ity of the wavefunctions that is imposed by the disorder 
present in the nanogram^ Random matrix theory of- 
fers a simple generic approach to modeling this fluctu- 
ating renormalization of the one-electron levels. In the 
presence of single-electron energy level renormalization 
the tunneling current /(V), at bias voltage V, is deter- 
mined by a complex master equation that contains the 
kinetics of all the allowed electron tunneling transitions — 
including also the renormalized energy levels. Those are 
given by the following Coulomb blockade threshold equa- 
tions: ns(V) > E™ x + U (for events involving unoccu- 
pied levels in the grain) and /zd(V) < E™ x + U (for 
events involving occupied levels in the grain); here n = 
represents the ground-state energy levels ( E® Xg = E a ). 
Hence, each time a new tunneling channel is opened a 
peak in the differential conductance is expected (see be- 
low, Sec. Hvjl. 

In summary, we have so far introduced the mechanism 
that creates electronic excitations, recast these excita- 
tions as a hierarchy of generations, stated that these ex- 
citations renormalize the single-electron energy levels in 
a complex manner, and discussed the effect of the renor- 



malization on the tunneling spectroscopy of the grains. 
In the following section we present a mean-field the- 
ory that can be used to calculate the renormalized one- 
electron levels of the excited nanoparticle. 



III. MEAN-FIELD THEORY OF THE EXCITED 
QUASI-PARTICLE STATES 

A neutral metallic grain with no valence electrons in its 
ground-state (|G)) is the starting point. We then adopt a 
quasi-particle picture consisting of electronic levels \tjj a ) 
with energy E a , and take as \G) the lowest energy no- 
electrons configuration [see Fig. ^a)]. By implementing 
a tight-binding Hamiltonian, E a and \ip a ) are calculated 
taking into account important microscopic features of 
the nanograin: geometry, structural disorder, and surface 
chemistry>! In a multi-orbital (s, p, and d) representation, 
and neglecting spin-orbit coupling, the Hamiltonian is: 

TC = E i c\ a Ci a + Tja,ia' {c\ a ' c ja + c ] Q c ia') ■ 

(1) 

i (j) labels the lattice site Ri (Rj) of atom i (j), c\ a 
{ci a ) creates (destroys) an electron on site i, while a (a') 
indicates the s, p or d orbital, sf and Tj a ^ a i are the 
Slater-Koster (SK) on-site and hopping parameters! 15 ' 16 

The metallic nanoparticles that are probed by electron 
tunneling spectroscopy are passivated with a thin Al- 
oxide layer and buried between two massive electrodes. 
Hence, to extract data about the morphology of these 
nanograms is quite difficult. However, their size — several 
hundreds of atoms — makes it reasonable to model their 
structure as a crystaline metal core, and a disordered 
superficial (shell) region that corresponds to the metal- 
oxide interface. The disorder present in the surface 
arises from the chemistry of the oxide layer in combina- 
tion with surface reconstructions that may occur during 
growthiii The model nanoparticles considered here re- 
sult from truncating a fee lattice to a volume V in a disc 
or hemisphere geometry^ In this case, the coordination 
numbers of atoms at the surface of the nanoparticle dif- 
fer from the coordination number of the bulk fee lattice. 
This is used to establish a criterion that defines the sur- 
face of the particle, and distinguishes between core and 
shelli Once the surface atoms are determined we ran- 
domly choose 50% of those sites to represent O while the 
other sites correspond to Al. This randomness is the only 
source of surface disorder we consider. 

We adopt as SK parameters in Eq. for the atomic 
sites in the core of the nanograin those parameterized by 
PapaconstantopoulosJ^ The aluminum and oxygen sites 
in the metal-oxide interface (shell), on the other hand, 
have different SK parameters as a consequence of charge 
transfer from Al to O. To find the on-site energies for 
the charged oxygen and aluminum atoms in the oxidized 
shell we combine (i) the Mulliken-Wolfsberg-Helmholz 
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(MWH) molecular-orbital approach*^ with (ii) results 
of classical molecular dynamics simulations by Camp- 
bell and co-workers^ that show the metal-oxide interface 
at the nanometer scale is mainly constituted of interca- 
lated O" 1 / 2 and AI+ 1 / 2 . The de gree of charge transfer 
v determines the on-site energies via the MWH theory: 
The molecular orbital [p) energies of atom m are em- 
pirically parameterized as E™ — E™(y) — — (A"V 2 + 
B™v + C™)}22i where v is the excess valence charge: 
v = (— )l/2 for Al (O). Hence, the SK on-site ener- 
gies for the atomic sites in the oxide layer are given by 
= ~(^ (0) /4 + {-)B^/2 + C^ 0) - A), with 
A = 1.289.% an off-set energy such that Ef(0) = e s Al , 
and the parameters A™, B™, and C™ for m = Al, O 
those of Ref . Finally, we determine the nearest neigh- 
bor O-Al hopping SK parameters by assuming an average 
separation between oxygen and aluminum atoms in the 
oxide of <Iai~o = 1.8Aj21 applying Harrison's model to 
obtain the two-center transfer integrals 2 ^ and then trans- 
forming them to find the SK hopping parameters^ By 
taking do-o = 3.0A— as average separation between 
atoms the above procedure leads to the nearest neighbor 
0-0 hopping parameter. We keep the SK hopping pa- 
rameters for Al +1 / 2 -Al +1 / 2 and Al +1 / 2 -Al the same as 
those for Al-Al. 

Let us now assume that electronic excitations of the 
no electrons within the grain are present due to the 
mechanism discussed in Sec [HI These excitations can 
be identified according to the number of electron-hole 
pairs generated over the ground-state \G). A generic n- 
electron-hole (n-eh) pair excitation is given by: \X n ) = 

\k<k<---k<;k>k>---k>) = (n? =1 ci>)(ns, =1 c fc <)|G} J 

K l m 

where kf (k^) labels an empty (occupied) quasi-particle 
state in the ground-state, and cl (c/s) creates (destroys) 
an electron in state \ipk}- As mentioned above, the exci- 
tation changes the electronic charge density within the 
grain from po{R), in the ground-state, to px n {R) — 
Pg{R) - eELi K>(^)| 2 ~ E:=i \^k<(R)\ 2 ] where R 
indicates an atomic site inside the grain, and \ipk{R)\ 2 
is the wavefunction amplitude of state \ipk) at site R. 
Therefore there is an induced {bare) charge density in 

the grain given by: Spx„(R) = Px n (R) — Pg{R) which 
is screened by the electrons present in the metallic 
grain. This screening can be modeled at different lev- 
els of complexity^* Here we implement for simplicity 
a static screening represented by an effective Thomas- 
Fermi (TF) dielectric constant e^ F = 1 + (qrF/q) 2 

(q TF = 4(37r 5 a 3 3 A/")~ 1/3 and q = 2/aAi with aAi = 
0.405A and Af = 0.18lA" 3 the bulk Al lattice param- 
eter and electronic density, respectively; as — 0.529A is 
the atomic Bohr radius). 25 Within this approach the net 
(screened) induced charge is: Spx„ (R)/e^ F . The choice 
of the the wave vector magnitude q that enters the TF 
screening is due to the typical length scale on which the 
electron wavefunction changes within the grain£ 

To investigate the effect of the electronic excitations 



{|A n )} in the energy spectra of the grain, we extend our 
tight-binding model by setting the on-site energy of or- 
bital a in atomic site Rj to: elf™ = Ej a + £j a . Here, 

is the renormalization energy introduced by the Coulomb 
interaction due to the electronic excitation. Within 
the empirical WHM approach discussed above^ Uj a = 
(2A J a qj + B^)e8px n {Rj)/^ rF \ an d accounts for the ex- 
cess charge [—e6px n (Rj)/e% F ] present in atomic site Rj 
due to the redistribution of charge inside the nanograin; 
qj is the ground-state charge present in site j. The 
second term on the right hand side of Eq. (J2J is the 
off-site Hartree contribution^ It should be noted that 
in the present model we do not consider any exchange 
effect — fine structure — that may distinguish between dif- 
ferent electron-hole spin configurations. Furthermore, 
we do not include scattering among different excitations 
{|X n )}. This approximation should be valid for the low- 
lying electronic excitations. 27 



IV. CASE STUDY: ALUMINUM NANOGRAINS 

We now address the effects of electronic excitations 
on the energy spectra of metallic nanograins by pre- 
senting detailed results for the response of the single- 
electron energy levels to the creation of the lowest-lying 
excitations that are accessible at low applied bias volt- 
age in a two-terminal device containing a disc-shaped 
Al nanograin coated with Al oxide. The grain's vol- 
ume V = 13.4 nm 3 , the drain-source capacitance ratio 
Cd/Cs = L6, and the charging energy U = iS; 28 here 
S = (4_E^ Z /3A/")V _1 = 6.3 meV is the average energy 
level spacing — predicted by a particle-in-a-box model of 
the grain — around the Fermi energy of bulk Al [Ep\. 
We consider three representative grains — A, B, and C — 
that differ in the specific realization of the disorder in 
their oxide coats. 



A. Ground state and excited energy levels 

The ground-state (\G)) single-electron energy struc- 
ture of nanograin A is shown in Fig. Efa); around its 
Fermi energy (Ep), and for a particular realization of 
disorder^ Figure Eta) also shows the renormalized en- 
ergy spectra for generation 1 that arise by reaching the 
threshold voltage Vg\ rd — V* h to inject an electron into 
level \ipE F +i) = the spectra are labeled by the par- 
ticular lch electronic excitation that is involved. For 
this device at the threshold bias voltage, the resulting 
singly charged nanoparticle can decay by emitting an 
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electron from any of seven different single-electron spin- 
degenerate orbitals into the drain contact. As previ- 
ously discussed, the different excitations renormalize the 
one-electron energy levels differently; and fluctuations in 
the magnitude of the energy renormalization, for a given 
level, are visible, ft should be noted that the renormal- 
ization of the energy levels is smaller than the average 
energy spacing of the ground-state single-electron energy 
levels^ (6) ~ 6.9 meV. Furthermore, for each excited 
state the average spacing between one-electron energy 
levels differs from (6) by less than 2%. Increasing the 
bias voltage to V^" 1 , at which it becomes possible for an 
electron to tunnel into level \"4>e f +2) = |2), leads to new 
leh-pair excitations. FigureEJb) shows in detail — within 
an energy interval of 36/5 — the excited single-electron 
levels around Ep + \ and Ep+2 that result from genera- 
tion 1. For comparison, the box in the lower panel of Fig. 
Hfb) shows the renormalizcd energy levels for generation 
2; only for Ep + i, and at V* . Clearly, the fluctuations in 
the renormalization of the energy levels depend strongly 
on the details of the excitation(s) that are involved. 

As mentioned above, disorder imparts a stochastic na- 
ture to the confined electronic wavefunctions, thus dra- 
matically affecting the renormalization fluctuations. The 
latter are visible in Figure |3 where the ground-state and 
excited electronic energy levels for different disorder re- 
alizations are shown: While for grain A the spectra show 
one of the renormalized levels of Ep+i and Ep+2 visibly 
separated from the others, for all these leh-pair excita- 
tions, this is not the case for B and C; see Fig. |3] In par- 
ticular, for grain C, the renormalized energy levels corre- 
sponding to Ep + 2, at V 2 th , superpose with the renormal- 
ized levels for Ep+3- It should be noted that the number 
of excitations is similar for the different nanograins. This 
results from the combined effect of (i) adopting the same 
drain-source capacitance ratio (Cd/Cs = 1.6) and charg- 
ing energy (U — 46) for all the grains, and (ii) having the 
volume (V) fixed while changing the disorder from grain 
to grain. The latter gives nearly the same average num- 
ber of single-electron energy levels in an energy interval 
of several 5 for each graini In particular, the number 
of available orbitals for generating excitations, at V-f 11 , 
is roughly: mt{[U + E F+1 - E F ][1 + {C D /C s )- 1 ]/6} = 
mi[6.5+1.625 (Ep+x — Ef)/6], where int(x) is the integer 
part of x. 

In conclusion, the response of the electronic structure 
to the excitations exhibits strong dependence on the en- 
ergy levels involved in the excitation, and disorder real- 
ization. 



B. Differential conductance 

Having calculated the renormalized energy levels let us 
turn now to the implications for the differential conduc- 
tance. We concentrate here on finding the energies — or, 
equivalently, bias voltage — at which tunneling resonances 
should be present rather than performing an actual cal- 



culation of the dl/dV spectrum. The latter is beyond the 
scope of this work as it would require to solve a compli- 
cate master equation describing the kinetics of the tun- 
neling transitions far from equilibrium. 

In general, resonances in dl/dV arise when new chan- 
nels for tunneling open as the bias voltage is swept. In the 
Coulomb blockade regime, three kinds of tunneling reso- 
nances are expected in the spectra of metal nanograins: 
Direct, charge- fluctuation^ and non-equilibrium^ The 
first two are present regardless of the ratio between 
the electronic relaxation and bottleneck tunneling rate, 
therefore they are regarded as equilibrium resonances. 
Direct resonances generally correspond to tunneling tran- 
sitions directly into (out of) an electronic energy level 
in the ground state of the nanograin. The charge- 
fluctuation resonances that are present at low bias 
voltage arise from having a finite probability for the 
nanograin to be (on average) negatively charged — excess 
of one electron £ These resonances are also non-trivially 
affected by non-equlibrium effects (see below), however, 
they were not included in the original AWARTi model. 
Here, we generically label them Q + . 

Figure 21 summarizes the bias voltages (V), presented 
in units of energy after converting to eV, at which tun- 
neling resonances should be present in the dl /dV spec- 
tra of grains A, B, and C; at low bias voltage. Reso- 
nances appear in groups: three in A, and two in B and 
C; only G#0 (ground-state) and G#l are shown. In 
these grains the Fermi level is doubly occupied, as men- 
tioned in Sec. ^ furthermore, the adopted capacitance 
ratio Cd/Cs is such that in these devices the onset of 
tunneling takes place when an electron is injected from 
the source electrode into the grain. Hence, the first res- 
onance in dl / dV appears when the bias voltage reaches 
V^ h . Upon increasing V, an additional non- equilibrium 
tunneling resonance appears in the spectrum each time 
a renormalizcd energy value of -Ef+i is reached. Further 
increase in the bias voltage results in new equilibrium 
and non-equilibrium resonances that arise from tunnel- 
ing into Ep+2 (at V^ 1 ) and its renormalized energy lev- 
els, respectively. Non-equilibrium resonances appear as 
satellites of the main, equilibrium, resonances due to tun- 
neling into Ep+i and Ep+2- This leads to group 1 and 
3 in A] 1 and, partially, group 2 in B; and group 1 and 
2 in C. It should be noted that the satellites involving 
level Ep+2, however, do not systematically appear above 
the equilibrium resonance. The latter would make the 
experimental identification of the equilibrium resonance 
difficult. Moreover, the number of such satellite reso- 
nances is bigger than in the neighborhood of V-f* 1 as the 
accessible leh-pair excitations at V^ h involve both levels 
Ep+i and Ep+2- Finally, at bias voltages greater that 
Vi h the finite probability to have the nanoparticle neg- 
atively charged leads to charge-fluctuation resonances® 
(equilibrium and non-equilibrium satellites) that are re- 
sponsible for group 2 (Q + ) in A, and part of group 2 in 
B. However, no Q + resonances are present in grain C as 
a consequence of its electronic structure below the Fermi 
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energy^i It is also noticeable that in grain C the clusters 
of resonances are less dense than in A and B although 
the number of leh-pair excitations is nearly the same in 
all grains, as discussed above. This shows the high sen- 
sitivity of the tunneling spectrum to the response of the 
electronic energy level structure to eh excitations, and to 
the disorder present in the grain. 

If we were to include in the above discussion of the 
dl/dV spectrum the renormalized energy levels corre- 
sponding to G#2 [see Fig. |2£b)] this would dramatically 
increase the number of non-equilibrium resonances. In 
fact, if the kinetics of the tunneling transitions is such as 
to admit a significant population of excited states then it 
should be possible to achieve much higher spectral densi- 
ties of tunneling resonances than observed to date, sim- 
ply by increasing the thickness of the tunnel barriers and 
hence the tunneling time relative to the relaxation time 
of the eh pair excitations. This kinetics is determined by 
microscopic parameters of the grain, and the nature of 
the electron-phonon scattering within it; which strongly 
depend on the volume of the grain (see Ref. l32|) and dis- 
order. Hence, there is no a priori (during the growth) 
control on the expected number of non-equilibrium reso- 
nances during charge transport in a device. Nonetheless, 
tunnelin g ex periments in gated devices (as those reported 
in Ref. l33|) . where the charging energy of the device 
may be changed with the applied gate voltage, should 
at least be able to address the increase of the number 
of non-equilibrium resonances as a function of the avail- 
able orbitals (N; see Sec. |Hj| to generate the electronic 
excitations. Such observations may be contrasted with 
the predicted number of excitations in generation n as a 
function of N; see Sec. |H] However, the details of the 
renormalization of the probed energy level may hinder 
such a comparison, [see Fig. |21and0] grain C, where the 
number of non-equilibrium resonances in the first cluster 
is smaller (nearly 1/2) than the number of accessible ex- 
citations, as a result of having some renormalized levels 
(spectrally) below Ep+i- The latter implies that these 
excitations contribute to determine the tunneling current 
at Vi h , rather than adding satellite resonances into the 
dl/dV spectrum.] 

To summarize, the high sensitivity of the ground-state 
and excited electronic energy spectrum to disorder is non- 
trivially inherited by the differential conductance spec- 
trum: (i) The number of neutral tunneling resonances 
(equilibrium and non-equilibrium) strongly depend on 
the response of the electronic energy levels to the ex- 
citations, (ii) Charge-fluctuation resonances appear in 
clusters that, depending on disorder, can overlap with 
the neutral non-equilibrium resonances, as shown in Fig. 
21 (grain B). This would increase the complexity of the 
tunneling spectrum and make the interpretation of ex- 
periments difficult. 

V. SUMMARY 

In summary, we have presented a systematic study 
of the tunneling-induced electronic excitations in non- 



magnetic metallic nanograins: (i) We have discussed the 
mechanism that creates electronic excitations and classi- 
fied them according to the resulting number of electron- 
hole pairs, which enabled us to introduce the concept 
of generations of excitations: A subspace of excitations 
with a specific number of electron-hole pairs. The gen- 
erations form a hierarchy, and the number of elements 
in a generation is determined by the applied bias voltage 
and characteristics of the device, (ii) We have proposed a 
general microscopic mean-field model that can be used to 
calculate the quasi-particle energy levels in the presence 
of a tunneling-induced excitation. Based on this model, 
we have presented detailed results for the response of 
the electronic structure of ultra-small aluminum grains 
coated with Al oxide to these excitations. However, the 
model should be applicable to most of the non-magnetic 
metallic nanograins currently probed in tunneling spec- 
troscopy experiments. Our results show that disorder 
present at the surfaces of the nanograins imparts a fluc- 
tuating character to the renormalization of the single- 
electron energy levels due to the excitations, and that 
this renormalization is smaller than the typical energy 
spacing between single-particle levels. Furthermore, the 
tunneling spectra of the nanograins consist of equilibrium 
and non-equilibrium resonances, which appear in clusters 
whose structure is dramatically affected by the high sen- 
sitivity of the ground-state and excited electronic energy 
structure to disorder. 

We have also shown that if the nonequilibrium reso- 
nances discussed here and in previous work^ are present 
in the tunneling spectra of metal nanoparticles at all, 
then their density in the spectra should vary greatly 
from sample to sample. This is because (a) the density 
of nonequilibrium resonances increases rapidly with the 
number of generations of excited states that are popu- 
lated in the nanograin and (b) this number of populated 
generations is in turn very sample dependent. 

Finally, we have suggested that tunneling spectroscopy 
experiments in gated ultrasmall nanograins may be useful 
to probe the variation in the number of non-equilibrium 
resonances as a function of the number of available or- 
bitals for tunneling. 
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FIG. 1: (a) Mechanism for creating electronic excitations, exemplified with the creation of a leh-pair excitation: from the 
source (S) electrode, an electron tunnels into energy level \ipF+i) = |1) of an originally neutral grain containing no conduction 
electrons in its ground state |G), with a charging energy U as represented in the plot; the grain gets negatively charged 
(no — > no + 1), and the single-electron energy levels are renormalized by an amount U; subsequently, an electron tunnels 
out — from, in this example, level \tpF-i) = |1) — to the drain (D) contact which leaves the nanoparticle in an electronic excited 
state: = c[ci\G). In the latter notation, c[ (cj) creates (destroys) an electron in state |1) (|1)). (b) G#l and G#2 are 

representative elements of generation 1, containing leh-pair excitations, and generation 2, corresponding to 2eh-pairs; details 
in text. 1 11; 13) = cjcjcjCjilG) is pictorially present in this example, (c) Schematic representation of the renormalization of the 
ground-state single-electron energy levels after generating excitation |1; 1) — see text. The renormalized energy levels are shown 
as dotted dashes. 



FIG. 2: (a) Energy levels for grain A (volume V = 13.4 nm 3 ) around the Fermi energy (Ef) of the ground-state grain in units 
of 5 — 6.3meV — the average electron-in-a-box energy level spacing — for the ground-state (|G), bold dashes) and different one 
electron-hole pair excitations (\Xi) = ctcv|G) = \k; k')\ k and k' indicate the (excited) electron and (excited) hole one-electron 
levels, respectively. Thin dashes.) that may be created after the applied bias voltage reaches the threshold for tunneling an 
electron into level |1) in the neutral grain from S: Vj* . (b) Details of the renormalization of levels \ipF+i) and \ipF+2) arising 
from the presence of leh-pair excitations at V* and V% — bias voltage to tunnel an electron into |i/>f+2) from S. The box 
marks the renormalized energy of single-electron level |1) that results from 2eh-pairs excitations. Full scale corresponds to 
energy interval of 35/5. 



FIG. 3: Ground state |G) (bold dashes) and excited single-electron energy levels (thin dashes) for different realizations of 
disorder — grains A, B, and C. Excited states belong to G#l (leh-pair excitations), at bias voltages Vl h and Vj' 1 (see text and 
caption of Fig. The renormalization of the energy levels shows a strong sample (disorder realization) variation. Notice the 
different energy scales in the plots. Ef is the Fermi energy of each grain, and 5 is the same as defined in Fig. |2] 



FIG. 4: Applied bias voltages (converted to energy: eV) derived from the calculated ground-state and excited energy spectra at 
which tunneling resonances would appear in the dl jdV spectrum of the studied devices (see text). Only generations (G#0) 
and 1 (G#l) are considered when finding the energies for direct and charge-fluctuation (equilibrium and non-equilibrium) 
resonances. Thin and bold (blue in on-line version) bars distinguish between energies associated with direct and charge- 
fluctuation resonances, respectively. Different heights mark equilibrium and non-equilibrium resonance energies: 3/2 and 3/4 
for the former, while 1 and 1/2 label the latter. In the top panel — grain A — the symbol Q + indicates a cluster of charge- 
fluctuation resonances (see text, and Ref. 8j). The * indicates nearly-degenerate values of eV. The ordinate scale is the same in 
all panels, however, notice the different range of bias voltages in each plot. U and 5 have the same meaning as in the previous 
figures. 
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